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Abstract 

In this paper, we present a general quantum computation compiler, which maps any given quantum algorithm 
CN to a quantum circuit consisting a sequential set of elementary quantum logic gates based on recursive cosine - 
O sine decomposition. The resulting quantum circuit diagram is provided by directly linking the package output 



(3jT) written in LaTex to Qcircuit.tex <http://www.cquic.org/Qcircuit>. We illustrate the use of the Qcompiler 
<^ package through various examples with full details of the derived quantum circuits. Besides its generality and 
simplicity, Qcompiler produces quantum circuits which reflect the symmetry of the systems under study. 
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Program summary 

Program Title: Qcompiler 

Journal Reference: 

Catalogue identifier: 

Licensing provisions: none 

Programming language: Fortran 

Computer: any computer with a Fortran compiler 

Operating system: Linux, Mac OS X 10.5 (and later) 

RAM: depend on the size of the unitary matrix to be decomposed 

Keywords: Quantum compiler, quantum circuit, unitary matrix, quantum gate, quantum algorithm 

Classification: 

External routines/libraries: Lapack 
Nature of problem: 

Decompose any given unitary operation into a quantum circuit with only elementary quantum logic gates. 
Solution method: 

This package decomposes an arbitrary unitary matrix, by applying the CSD algorithm recursively, into a series of 
block-diagonal matrices, which can then be readily associated with elementary quantum gates to form a quantum 
circuit. 

Restrictions: 

The only limitation is imposed by the available memory on the user's computer. 
Comments: 

This package is applicable for any arbitrary unitary matrices, both real and complex. If the unitary matrix is real, its 
corresponding quantum circuit is much simpler with only half number of quantum gates in comparison with complex 
matrices of the same size. 
Running time: 

Memory and CPU time requirements depend critically on the size of the unitary matrix to be decomposed. Most 
examples presented in this paper require a few minutes of CPU time on Intel Pentium Dual Core 2 Duo E2200 @ 
2.2GHz. 
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1. Introduction 



Quantum computing exploits the nature of the quantum world in a way that promises to solve problems 
which are intractable using conventional computers Q3-(3]|. At the heart of a quantum computer lies a set of 
qubits whose states are manipulated by a series of elementary quantum logic gates, namely a quantum circuit, 
to provide the ultimate computational results. In our earlier work fl4), we developed a highly efficient quantum 
computation simulator to assist on the analysis of complicated quantum circuits comprised of qubit and qudit 
quantum gates. In this paper, we present a quantum computation compiler which maps any given quantum 
algorithm to a quantum circuit consisting a set of elementary quantum logic gates. 

In the seminal papers by Barenco et al and Deutsch et al [ElO, it was proven that any arbitrarily complex 
unitary operation can be implemented by a quantum circuit involving only one- or two-qubit elementary 
quantum logic gates. Earlier studies applied the standard triangularization or QR-factorization scheme with 
Givens rotations and Gray codes to map a quantum algorithm to a series of elementary gate operations [S|5]- 
0. Subsequently, a more efficient and versatile scheme based on the cosine-sine decomposition was proposed 
and utilized ItSUTBl. More recently, De Vos et al |[T4~l[T5l examined another decomposition scheme, namely 
the Birkhoff decomposition, which was found to provide simpler quantum circuits for certain types of unitary 
matrices than the cosine-sine decomposition. However, the Birkhoff decomposition does not work for general 
unitary matrices. 

In this work we develop a general quantum compiler, named the Qcompiler, based on the cosine-sine 
decomposition scheme, because it works for arbitrary unitary matrices, it is extremely adaptable, and the 
core CSD code is now available as part of the LAPACK package. Furthermore, we have managed to half 
the number of quantum gates in the circuit, if the quantum algorithm involves only real unitary matri- 
ces. For the Qcompiler, the input file contains a user specified unitary matrix U. Qcompiler applies the 
cosine-sine decomposition (CSD) recursively and compiles U into a quantum circuit with a sequence of 
controlled/uncontrolled rotation and phase gates. Its output contains the complete information on these el- 
ementary quantum gates, including a separate LaTex document which can be directly linked to Qcircuit.tex 



<http://www.cquic.org/Qcircuit> to produce the resulting quantum circuit diagram. 

This paper is structured as follows. Section 2 describes the recursive cosine-sine decomposition scheme, 

which maps an arbitrary unitary matrix U into a quantum circuit with only one- and two-qubit logic gates. 

In Section 3, we provide a more cost-efficient method for real unitary matrices, which significantly reduces 

the number of quantum gates in the final circuit. In Section 4, we discuss the general structure as well as the 

usage of the Qcompiler package. In Section 5, we present various examples with full details and discussions 
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of the resulting quantum circuits. Section 6 concludes the paper. 



2. The recursive CSD scheme: General Unitary Matrices 

The Cosine Sine Decomposition (CSD) algorithm as described in [fT6lfT7l decomposes an arbitrary 2" x 2" 
unitary matrix U as the following 



where 



U 



( \ t \( \ 
u 



v v > 



C S 

-s c 



x 

I y. 



(i) 



C = diag (cos #/) = 

/=1,...,2"-' 



cos 6\ 

COS 6 2 



COS 9 2 n-i 



and 



S = diag (sin#/) 

1=1,.. .,2"-' 



sin#! 



sin 09 



sin Q 2 n-\ 



are 2" 1 x 2" 1 diagonal matrices, while u,v, x and y are 2" 1 x 2" 1 unitary matrices which can be further 
decomposed by the CSD algorithm recursively into a string of block diagonal unitary matrices 0. 

To show the general structure of the decomposition, we denote all sub-matrices as u',. ignoring their 
explicitly different numerical values, where i indicates the level of recursion and dim the matrix dimension. 

At the first level, we have 



U 



2"-'x2" 



2"-'x2" 



2"^ 1 x2"- 1 / 



-1 

I 

-s 1 



c 1 s 

2"- 1 x2"- 1 2"-'x2"- 1 



c 1 

V "2"-'x2"- 1 2"- 1 x2"- 1 7 



(2) 



= U l A l U l . 



At the second level, 



M 2"- 1 x2"- 1 ~~ 



2"- 2 x2"- 2 



U 



2"~ 2 x2"- 2 J 



r 2 s 2 

2"- 2 x2"- 2 2"- 2 x2"- 2 

s 2 r 2 

V 2"- 2 x2"- 2 2"- 2 x2"- 2 



2«-2 x2 «-2 



u 



2 n-2 x2 «-2 



2"-'x2"- 1 



1 

2«-i X 2«-iy 



2"- 2 x2"- 2 



2"- 2 x2"- 2 



U 



2"- 2 x2"- 2 



U 



2«-2 x 2n-2/ 



r 2 s 2 

2"- 2 x2"- 2 2"- 2 x2"- 2 

s 2 r 2 

2 n - 2 x2"- 2 2"- 2 x2"- 2 



c 2 s 2 

<-2"- 2 x2"- 2 J 2'"- 2 x2"- 2 

-s" 2 r 2 

2"~ 2 x2"- 2 2"- 2 x2"" 2 ^ 



U 



2 n-2 x2 n-2 



l 2"- 2 x2"~ 2 



u 



2"- 2 x2"- 2 



U 



2 n - 2 x2 n - 2 > 



and 



= U 2 A 2 U 2 , 



U = U l A l U l = U 2 A 2 U 2 A l U 2 A 2 U 2 . 



At the ?'th level of recursion, the matrix U' is decomposed as the following 



i/ M = u'a'u' 



where 



2"-'x2"-' 



u l 



2"-'x2"- 



(6) 



2"-'x2"-'/ 



and 



C S' 

2"-'x2"-> 2 n -'x2"-' 

- S ' C' 

2 n -'x2 n -' 2"-'x2"-' 



c 1 s l 

2 n -<x2 n -' 2"-'x2 n -' 

- S ' C' 

2"~ i x2 n ~ i 2"-'x2"-'J 



At the end of the recursive process, i.e. the nth level, we obtain 



U =U l A l U { 



= U 2 A 2 U 2 A l U 2 A 2 U 2 



-- U 3 A 3 U 3 A 2 U 3 A 3 U 3 A l U 3 A 3 U 3 A 2 U 3 A 3 U 3 



(2"-\ 



n 

\ P =\ 



n a i(p) 



U n p A p 



rjn 



(7) 



(8) 



where p marks the position of the matrix sequence, i is implicitly determined by p, A is given by Eq. (|7j), and 



U n = diag (u n k ) = diag (expOV M ))- 

k=l,...,2" ' k=\,...,2" 



(9) 



The above described recursive CSD scheme works for 2" x 2 n unitary matrices. For an arbitrary NxN unitary 
matrix U, where 2 n ~ l < N < 2", we add an identity matrix to U to form a new unitary matrix 



W2"x2" 



NxN 



I(2"~N)x(2"-N) ) 



(10) 



and then apply the recursive CSD decomposition to W as described above. 

The decomposed matrices given by Eq. ([8]) can be directly related to elementary quantum gates, in par- 
ticular, the phase gate <D and the controlled rotation gate (31 19), where the rotation operation is defined as 



Ra(p) = exp(m • cr^) = / cos ^ + z'a • <x sin ^ . 



If the rotation axis is y or z, we have 



(ID 



and 



respectively. 



R y (p) = R y (29) = exp(icr y 9) = I cos 9 + icr y sin 9 



( \ 
cos 9 sin# 



- sin 9 cos 9 



(12) 



R z (p) = R z (2cp) = exp(/o%0) = / cos (p + icr z sin (p 

expO'0) 
exp(-/0)^ 



(13) 



To establish a mapping to elementary quantum gates, Mottonen et al [|9]| inserted identity matrices / = 
P l p P \P l p P) ^) after each A 1 ^ in Eq. (8 ) with which P'p P) commutes, and thus 



Un A i( P ) = un A m P m^pi^ = U n p P* p) A* p \pWf . 



(14) 



The final decomposition becomes 



/2"-l 



u = 



fjn 
u 2" ' 



(15) 



where U n p = (P i( p p _~Su n p , except U n x = U n v and B i( p p) = U n p P i( p p) is a 2" x 2" diagonal unitary matrix. With a set 
of specified by solving a set of linear equations, Bp P) will have the required symmetry to be equivalent to 



gate 



C n ~ l R z {i; !,...,«), 



(16) 



where i denotes the target qubit with n - 1 control qubits. Similarly, A 1 can be mapped to gate 



C n ~ l R y {i; l,...,z - 1,/+ l,...,n), 
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(17) 



and WL is equivalent to a series of gates 



(18) 



0(l;)<g)/ 2 *-i x2 »-i, 

R z (l;)®I 2 n-l x2 n-l, 
C l R z (2\l)®I 2 n-2 x2 ,,-2, 

C n - 2 R z (n-\;l,...,n-2)®I 2x2 , 
C n ~ l R z (n\\,...,n- 1). 



As an example, we decompose a 2 3 x 2 3 unitary matrix by applying the CSD scheme recursively, i.e. 

U = 5 3 A]5^2 5 3 A 3 5 l A l 5 3 A 3 5 2 A 2 5 3 A 3^3 ; (19) 



where 

A Y p (p = 4) = C 2 R y (\;2,3), 
A z p (p = 2,6) = C 2 R y (2; 1,3), 
A 3 (p= 1,3,5,7) = C 2 i? v (3; 1,2), 

(20) 

5;(^ = 4) = C 2 ^(1;2,3), 
B 2 p (p = 2,6) = C 2 R Z (2; 1,3), 
fij(p = 1,3,5,7) eeC 2 ^(3; 1,2), 

and t/g = <D(1; ),/? z (l; ), C l R z (2; 1), C 2 i? z (3; 1, 2). For simplicity, the identity matrices are omitted in the 
above expressions. The equivalent quantum circuit for U given by Eq.[19]is shown in Fig. [TJ 



For a 2 4 x 2 4 unitary matrix, we have 



U =5tA|^A^^ J B 2 A 2 J B 4 A 4 J B 3 A 3 J B 4 A 4 J B>Ai 

^ A 4 D 3 A 3 d 4 a 4 t}2 a 2 d 4 a 4 d 3 a 3 D 4 A 4 



(21) 



5 9 A 9 5 10 A 10 5 11 A 11 5 12 A 12 5 1 3A 13 5 14 A 14 ^ 15 A 15 ?7 16 . 
The corresponding quantum circuit is shown in Fig. [2} 



8 



- nu - 



- ul 



Ry 



t3 



R- 



Ry - R, 



Rr 



R- 



RyUR, 



-a- 



-a- 



R- 



Ul A] B> A\ B\ A\ B\ A\ B\ A\ B\ A\ B\ A\ B\ 



- Ul - = 



O - R 



R- 




R, - 



Figure 1: Quantum circuit for an 8-by-8 complex unitary matrix. 
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Figure 2: Quantum circuit for a 16-by-16 complex unitary matrix. 



3. The recursive CSD scheme: Real Unitary Matrices 



If the 2" x 2" unitary matrix U is real, we have u n pk = exp(i(p p k ) = 1 or -1 from Eq.|9|and Up becomes a 



diagonal matrix consisting of only 1 or -1 . In this case, we can insert the identity matrix / = U n p (U n p ) 
after each Ap P) in Eq. [8j i.e. 



u;u; 



(22) 



p 

where U n p = \\ U n q . The decomposition for a real unitary U then becomes 



U = 



(2 n -\ 

n 

\ P =i 



A'(P) 



fjn 



(23) 



where A 1 ^ = U n A l p P> U n . Note for general unitary matrices, it requires the solution of a set of linear equations 



1 i(p) 



to determine each inserted P'i p) in Eq. 14 [9J. Here, for real unitary matrices, we only need to calculate 



the product of U\ , U\, U" to obtain the inserted U" Furthermore, we end up with only half number of 



decomposed matrices in comparison with Eq. 15 Since U" is a diagonal matrix with only two possible 



values, 1 or -l,Ap can be readily mapped to the gate C n l Ry(i; 1, i— 1, i+ 1, n). The gates for A'p" and 



n(p) 



Ap P) are almost the same except for the signs of some rotation angles. 



The last matrix [/?„ can be written as 



U'l„ = M diag (D mik ) g) 7 2 m-i x2 m-i (24) 

A J- J — 1 Of7-m 



m=l 



/t=l,...,2" 



where D m>fc are either ( o ? ) or ( o -°i )' w ^ m me former being an identity gate and the later a IT gate. Therefore, 
U\ n is equivalent to a subset of the following quantum gates: 

n(l;)<8>/ 2 , 1 -i x2 «-i 

C'n(2; 1) <g> I 2 n-2 x2 n-2 



C n - 2 U(n- l;l,...,7i-2)®/ 2x2 
C n_1 n(n; 1). 



As an example we decompose a 16-by-16 real unitary matrix U and get 

iy _ 14 T3 r4 T2 r4 r3 T4 ri p 74 72 74 73 74 fj4 

Its corresponding circuit is, as shown in Fig.[3j much simpler than that for a 16-by-16 complex unitary matrix. 
The number of gates is reduced to half. 

10 




u 4 

u 16 



A 4 A 3 A 4 A 2 A 4 A 3 
A 15 A 14 A 13 A 12 A U A 10 



A a An 



A 1 A 6 A 5 



A 4 A 2 
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4. Overview of the Software 



Figure 3: The circuit of a 16-by-16 real U. 



The Qcompiler package, written in Fortran, consists of 1 main program and 12 subroutines. This 
package utilizes recursively Sutton's CSD subroutines IfPTTl . which is now part of Lapack available at 
http://www.netlib.org/lapack under a permissive free software license. 

As described in Section 3, the Qcompiler package decomposes U as 



U = 



( /2"-l 

P =\ 

2"-l 



i(p) 



jjn 



if U is real, 



p | ^ 2" 

I! B^A^) U" if U is complex. 



(26) 



It then generates an output file containing a complete description of each quantum gate in the circuit. The 
principal flow chart of Qcompiler is shown in Fig. [4] with all subroutines calls summarized in Fig. [5} 
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Real Case 
Main Program 
CYGR_CSD 

i 



Read in the real 
unitary matrix 
CYGR_READF 
i 




Create a useful 
table called "index" 
CYGJNDEXTABLE 



CSD recursively 
CYGR_BLKCSD 



Start 

Main Program 
CYG CSD 




1 



Complex Case 
Main Program 
CYGC_CSD 

i 



Read in the complex 
unitary matrix 
CYGC_READF 



Create a useful 
table called "/ndex" 
CYGJNDEXTABLE 



CSD recursively 
CYGC_BLKCSD 



Process the CSD result in the real 
case to match each decomposed 
matrix with a quantum gate 
CYGR_CUTGATE 

T 



IE 



Process the CSD result in the complex 
case to match each decomposed 
matrix with a quantum gate 
CYGC_CSDPHASE 
T 



Create the output file (details 
of each gate in the circuit) 
CYGR_BLKCSD 




Create the output file (details 
of each gate in the circuit) 
CYGC_BLKCSD 







End 



Figure 4: Flow chart for the Qcompiler package. 



CYGR CSD 



CYGR READF 



CYGR BLKCSD 



CYGR CUTGATE 



CYGR WRITEF 



CYG_CSD 



CYGJNDEXTABLE 



CYGC CSD 



CYGC READF 



CYGC BLKCSD 



CYGC CSDGATEU 



CYGC CSDPHASE 



CYGC WRITEF 



CYGC_COEFF 



Figure 5: Flow chart for all subroutine calls. 
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5. Example 



5.1. Complex unitary matrix 

An 8-by-8 complex unitary matrix U is generated randomly as the following: 



U 



0.6501-0.3423i 
-0.2339-0.1113; 
-0.1885-0.1242; 
-0.2882+0.1223; 

0.0801-0.1293i 
-0.1002+0.0211; 
-0.2579-0.1785; 
-0.3122+0.1212; 



-0.0792-0. 1132i 
0.5429+0.0637; 
-0.1468-0.2134; 
-0.2745+0.1175; 
-0.4807+0.1042; 
-0.1428-0.4283; 
-0.0825-0.2606; 
0.0733+0.0072; 



-0.1411-0.0356; 

-0.0551-0.0559; 

0.6018-0.1214; 

-0.0532-0.2853; 

-0.3721+0.0264; 

-0.2636+0.1747; 

-0.2869+0.1512; 

-0.4026+0.0719; 



-0.1727+0.1594; 
-0.0454-0.1349; 
-0.5436-0.1291; 
0.2171-0.4183; 
-0.2358-0.3746; 
-0.1363+0.2767; 
-0.1106-0.2008; 
0.0101+0.2156; 



-0.2945-0.2656; 
-0.0212+0.0956; 
-0.1374+0.0600; 
-0.5067-0.4492; 
0.4089+0.0434;' 
-0.3588-0.0596; 
0.1490+0.0484;' 
-0.0048-0.162l! 



-0.3970+0.0475; 
-0.1494-0.0384; 
-0.1848-0.1079; 
-0.0320-0.0010; 
0.1480+0.1193; 
0.5210-0.2426; 
-0.3993-0.1107; 
-0.4437-0.1859; 



-0.1703-0.0385! 

-0.5056-0.1816; 

0.0026+0.1374; 

-0.0223+0.0879; 

-0.2750+0.1932; 

-0.0680-0.2148; 

0.5480-0.3288; 

-0.1984+0.2056; 



-0.0757+0.1065; \ 

-0.5333-0.0362; 

0.0116-0.3250; 

-0.1129+0.1611; 

-0.1832-0.2196; 

-0.1792-0.2124! 

-0.1874+0.1716; 

0.4193-0.3917; 



(27) 



Recursive CSD decomposition gives U = B^^B^A^B^^B^^B^A^B^A^B^A^U^. Their corresponding quan- 
tum gates are detailed in the output file produced by Qcompiler, as explained below. 



Output of the package 
GATEPHASE 

l; 

0.1586 



GATEZ 

l; 

-0.9734 



GATEZ 
2;1 

-0.0856 0.4745 



Matrix 



GATEZ 
3; 1,2 

0.0975 -0.0388 -0.4445 0.0461 



0(1;) 

O = 2tt- 0.1586 



Gate 



2(f>' j = 2n ■ (-0.9734) 



-R, - 



C l R z (2; 1) 

2<p\ = 2n ■ (-0.0856) 
2<p\ = 2n • 0.4745 



- R 7 -=- Rl-R 2 - 



C 2 R Z (3;1,2) 
2<p\ = 2n- 0.0975 
2<p' 2 = 2n ■ (-0.0388) 
2^3 = 2n ■ (-0.4445) 
20' 4 = 2n- 0.0461 



R z - 



- R - R - Rl - Rz - 
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Output of the package 



Matrix 



Gate 



GATEY 
3; 1,2 

0.0758 0.2318 0.4094 0.2817 



C%(3;1,2) 
2#i = 2n • 0.0758 
20 2 = 2n- 0.2318 
2(9 3 = 2tt • 0.4094 
20 4 = 2n- 0.2817 



: 3 



*> 2 



3 



GATEZ 
3; 1,2 

-0.5384 0.1700 -0.6665 -0.2091 



C 2 i? z (3;l,2) 

20! = In • (-0.5384) 

20 2 = 2tt- 0.1700 

203 = In • (-0.6665) 

204 = 2n -(-0.2091) 



-Rl-Ri-Rl-Rt- 



GATEY 
2; 1,3 

0.1573 0.3405 0.1019 0.4831 



C 2 i?,(2;l,3) 
20! = In- 0.1573 

20 2 = In • 0.3405 

203 = In- 0.1019 

20 4 = 2tt- 0.4831 



Ry--- 



TTTT 

pi p2 p3 p4 
— K y - K y - K y " K y ~ 

JTTT 



GATEZ 
2; 1,3 

-0.0876 -0.2271 0.2202 0.2315 



C 2 R z (2;l,3) 

201 = 2n • (-0.0876) 

202 = 2n -(-0.2271) 

203 = 2n • 0.2202 

204 = 2tt- 0.2315 



R- 



Rl 



Ri 



Rt 



Rt 



-X JTIJ 



GATEY 

3; 1,2 

0.1231 0.3567 0.2444 0.2683 



C 2 i?,(3;l,2) 

201 = 2n- 0.1231 

20 2 = 2tt • 0.3567 

20 3 = 2n • 0.2444 

20 4 = 2n • 0.2683 



: 3 
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Output of the package 



Matrix 



Gate 



GATEZ 
3; 1,2 

-0.3146 0.4879 



0.5310-0.1047 



C 2 i? z (3;l,2) 

201 = 2n -(-0.3146) 

202 = 2tt • 0.4879 

203 = 2n -(-0.5310) 

204 = 2tt -(-0.1047) 



~1 



- Rl - R 2 - Rl - Rt - 



GATEY 
1; 2, 3 

0.0901 0.1279 0.2752 0.3789 



A 1 



C 2 R y (l;2,3) 

26 1 = 2n • 0.0901 

26 2 = 2n- 0.1279 
20 3 = 2tt • 0.2752 
26> 4 = 2n • 0.3789 




GATEZ 

1;2,3 

0.5219 0.5019 0.2641 -0.0838 



C 2 i? z (l;2,3) 
20! = 2n • 0.5219 

20 2 = 2n ■ 0.5019 

203 = 2n ■ 0.2641 

204 = 2n • (-0.0838) 



R- 



Rl 



Rl 



Rl 



Rt 



GATEY 

3; 1,2 

0.1316 0.2352 0.2069 0.3275 



C 2 ^(3;l,2) 
2d x = 2n- 0.1316 
26 2 = 2n • 0.2352 
26> 3 = 2n ■ 0.2069 
26> 4 = 2n • 0.3275 



"2 



GATEZ 
3; 1,2 
-0.5367 



0.3010 -0.0056 0.5324 



C 2 i? z (3;l,2) 

20! = 2n • (-0.5367) 

202 = 2n • (-0.3010) 

203 = 2n • (-0.0056) 

204 = 2n • 0.5324 



~1 



R- 



Rl 



Rl 



Rl 



Rt 



15 



Output of the package 



Matrix 



Gate 



GATEY 
2; 1,3 

0.1898 0.4103 0.2375 0.3869 



C 2 R y (2;l,3) 
29 i = 2n- 0.1898 

20 2 = In • 0.4103 

20 3 = In • 0.2375 

20 4 = 2n ■ 0.3869 



-I 



=- Rl - Ri - Rl - Rt - 



GATEZ 
2; 1,3 

-0.1607 -0.5540 0.0320 0.8280 



C 2 R Z (2;\,3) 
20j =2n -(-0.1607) 
20 2 = 2tt • (-0.5540) 
2^3 = 2tt • 0.0320 
204 = 2n ■ 0.8280 



=-l # H H H /c 



GATEY 
3; 1,2 

0.1077 0.4126 0.3956 0.2628 



C 2 i?,(3;l,2) 
20! = 2tt- 0.1077 

20 2 = 2n • 0.4126 

20 3 = 2n • 0.3956 

20 4 = 2n ■ 0.2628 




GATEZ 
3; 1,2 

-0.0073 -0.6907 -0.3105 -0.4417 



C 2 R Z (3;1,2) 

20! = 2tt • (-0.0073) 

202 = 2n ■ (-0.6907) 

203 = 2^ -(-0.3105) 

204 = 2n • (-0.4417) 







3 C 

) < 


) 1 


i — i 


1 — 

( — 


- 


% - - R 








4 

z 



The final quantum circuit of the complex unitary matrix U given by Eq. [27] is shown in Fig. [6] 



- «W - 



- <D - R 



R- 



-a- 



Ry - R, 



-a- 



-a- 



7?. 



t2 



# y - R, 



^7 *7 ^6 *6 ^5 5 5 ^4 *1 ^3 *3 ^ 



-a- 



-a- 



^ ^1 



#z 

5? 



Figure 6: The circuit of U, the random 8-by-8 complex unitary matrix. 
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It requires 8 gates for 4 x 7 R y gates for A' p , and 4x7 R z gates for 5^, thus 64 gates in total to build the 
equivalent quantum circuit. 

5.2. Real unitary matrix 

An 8-by-8 real unitary matrix U is generated randomly as the following: 



U = 



-0.2991 
-0.3284 
-0.3766 
-0.3183 
-0.3415 
-0.4157 
-0.3908 
-0.3427 



0.1387 

-0.4431 

-0.7104 

0.1523 

0.0206 

0.4240 

0.1596 

0.2257 



0.5667 
-0.0415 

0.1156 
-0.7265 
-0.0599 

0.2668 
-0.2301 

0.0914 



0.0990 
-0.1607 
-0.1142 
-0.1032 
0.9064 
-0.2504 
-0.0746 
-0.2255 



0.5322 
-0.4763 
0.2958 
0.5063 
-0.1500 
-0.1871 
-0.2757 
-0.1124 



0.1225 
-0.0154 

0.0275 
-0.0310 
-0.1088 

0.2838 

0.4422 
-0.8337 



0.4933 
0.3696 
-0.3785 
0.0024 
-0.1481 
-0.5707 
0.3534 
0.0655 



-0.1371 \ 

-0.5519 

0.3094 

-0.2826 

-0.0394 

-0.2689 

0.6057 

0.2456 



(28) 



Recursive CSD decomposition gives U = A^A^A^A^A^A^A^U^. Their corresponding quantum gates are 
detailed in the output file produced by Qcompiler, as explained below. 



Output of the package 



Matrix 



Gate 



GATEPI 

l; 

Y 



GATEPI 
2; 1 
YN 



GATEPI 
3; 1,2 
Y YN Y 



GATEY 
3; 1,2 

-0.3188 0.4501 



-0.4725 0.2393 



n(i;) 



c l n(2- 1 ) 



C 2 n(3;(l,2) 00 ) 
C 2 n(3;(l,2) 01 ) 
C 2 n(3;(l,2) n ) 



-Ti- 



ll 



7T 



Note: This part specifies the subset of gates: 



n 




n 



n - 



C 2 R y (3;\,2) 

29 l = 2n- (-0.3188) 

26 2 = 2n ■ 0.4501 

26 3 = 2n ■ (-0.4725) 

26 4 = 2n ■ 0.2393 



3 



Ry 



R} 
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Output of the package 



Matrix 



Gate 



GATEY 
2; 1,3 

0.0624 0.4624 0.0128 0.4629 



A 2 



C 2 i?,(2;l,3) 

201 = 2n ■ 0.0624 

20 2 = 2n • 0.4624 

20 3 = 2n ■ 0.0128 

20 4 = 2tt • 0.4629 



=- Rl -Rt -Rl -Rt - 



X- JTTI 



GATEY 
3; 1,2 

-0.1026 0.4866 -0.0855 -0.2735 



C%(3;\,2) 

20! = In -(-0.1026) 

20 2 = 2n • 0.4866 

20 3 = 2n • (-0.0855) 

20 4 = 2n • (-0.2735) 



Rr 



Rt 



GATEY 

1;2,3 

-0.0152 0.1125 -0.3422 -0.4830 



A 1 



C 2 i?,(l;2,3) 

20! = 2n • (-0.0152) 

20 2 = 2tt- 0.1125 

20 3 = 2n • (-0.3422) 

20 4 = 2n • (-0.4830) 



- R y - - Rl - Rl - Rl - Rt - 



GATEY 
3; 1,2 

0.3353 -0.3157 -0.2717 0.1060 



C%(3;\,2) 

201 = 2n • 0.3353 

20 2 = 2tt -(-0.3157) 

20 3 = 2n -(-0.2717) 

20 4 = 2n- 0.1060 



RyU 



- Rl -Rt - Rl - Rt - 



GATEY 
2; 1,3 

0.2953 0.3166 0.2673 0.4561 



A\ 



C%(2;\,3) 
20! = 2n ■ 0.2953 

20 2 = 2n- 0.3166 

20 3 = 2n ■ 0.2673 

20 4 = 2n- 0.4561 



Ryl 



=- Rl - Rt - Rl - Rt - 



-J 
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Output of the package 




Gate 



GATEY 
3; 1,2 

0.2999 -0.2617 -0.4684 -0.281 1 



C%(3;\,2) 
20! = In ■ 0.2999 

20 2 = 2n • (-0.2617) 

20 3 = In • (-0.4684) 

20 4 = 2n- (-0.2811) 



*> 2 



Rl 



R, 



The final quantum circuit of the real unitary matrix U given by Eq. 28 is shown in Fig. |7} which is signifi- 
cantly simpler than that shown in Fig.[6]for a complex unitary matrix of the same size. 



n 



n | — o — »i — »>- 



ft 



n - n - n - R y 



Rr 



R, 



T "^f - f 



Ry 



Ry 



Ry- 



X3 VI X3 XI 43 42 A"3 
A n /i g /i 5 /i 4 /i 3 /i 2 /lj 



Figure 7: The circuit of f/, the random 8-by-8 real unitary matrix. 



It requires 5 IT gates for Ul, and 4x7 R y gates for A' thus 33 gates in total to build the equivalent quantum 



circuit, significantly simpler than its complex counterpart. 



5.3. Random walk on a square graph 

Quantum walks have recently been explored for their non-intuitive dynamics, which may hold the key to 
radically new quantum algorithms lfT8l - l2"0ll . Douglas and Wang 11211 recently presented efficient quantum 
circuits for several families of highly symmetrical graphs. Here we use Qcompiler to build quantum circuits 
for quantum walk on general graphs. 

The walker's quantum state is given : 



N 



i=l jeS 



(29) 



where /V is the number of nodes of the given graph, \i) represents the ith node state, \j) the coin state which 

connects the ith and jth node, and S is determined by the graph adjacent matrix. 
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To realize one step quantum walk, we first apply a coin operator C = diag (Q) on the walker, where 

i=l,...,N 



Ci = 2^2^c) k \Uk){iJ\. 

jeS keS 



(30) 



In this work, we choose the Grover Coin ll22Tl so that the quantum walk can be represented by a real matrix. 
We then apply the translation operator 



f = J]J]\j,i)(iJ\. 



i=\ jeS 



(31) 



The unitary matrix for the complete step is U = TC. 

For a simple square graph, shown in Fig. [8j the quantum walk operator is 



, 1 

i 



U = TC = 



(0 1 
1 



1 

1 



1 

1 



1 
1 



/OOOlOOOOx 
1 
10 
10 
10 
10 
10 

^00001000' 



Figure 8: The square graph 

The output file of Qcompiler is given as the following with its quantum circuit shown in Fig.|9j 



Output of the package 


Matrix 


Gate 


GATEPI 

l; 

Y 




n(i;) 


-& 


GATEPI 
2; 1 
YN 




C l U(2; 1 ) 


- 7T - 
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Output of the package 



Matrix 



Gate 



GATEPI 
3; 1,2 
N Y YN 



C 2 n(3;(l,2) ()1 ) 
C 2 n(3;(l,2) 10 ) 



7T 



7r 



Note: This part specifies the subset of gates: 



7T 




7T 



7T 



GATEY 
3; 1,2 

0.5000 0.0000 0.0000 -0.5000 



C 2 i?,(3;l,2) 

20i = 2n • 0.5 = n 

26 4 = 2n ■ (-0.5) = -7T 



~1 



- fl! - - 



GATEY 
2; 1,3 

0.0000 0.5000 0.0000 0.5000 



A 2 



C 2 i? y (2;l,3) 

26 2 = 2n • 0.5 = 7T 

2fl 4 = 2tT • 0.5 = 7T 



*> 4 



GATEY 
3; 1,2 

-0.5000 -0.5000 0.0000 0.0000 



C 2 i?,(3;l,2) 

20 x = 2n • (-0.5) = 

20 2 = 2n ■ (-0.5) = 



-7T 

—n 



-2 



GATEY 
1; 2, 3 

0.0000 0.0000 -0.5000 0.5000 



4 1 

A 4 



C%(\;2,3) 

20 3 = 2n • (-0.5) = -n 

20 4 = 2n • 0.5 = 7T 



R} 



GATEY 
3; 1,2 

0.0000 0.5000 0.5000 0.0000 



C%(3;1,2) 

20 2 = 2n • 0.5 = n 

20 3 = 2n • 0.5 = n 



r 



Rl 



Rl 
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Output of the package 



Matrix 



Gate 



GATEY 
2; 1,3 

0.0000 -0.5000 0.0000 0.5000 




C 2 tf v (2;l,3) 



20 2 

20 4 



In • (-0.5) = -n "J 
2n • 0.5 



n 





— • 


- R2 y 




ii — 





GATEY 
3; 1,2 

0.0000 -0.5000 -0.5000 0.0000 



C 2 R y (3-\,2) 

26 2 = In ■ (-0.5) = -n 

26 3 = 2n ■ (-0.5) = -n 



- R y H R, - 



r 



n — o- 



n 



n 



n 



Rt 



4 2 

A 6 



Rl 



Rl 



— o ii ii ii it 



A 1 



3 



^3 



Rz 



R 



Rl 



Figure 9: The circuit of U, the random walk evolution on the square graph (Fig.[8]l 



Compare to the random 8-by-8 real matrix given by Eq. 28 , which requires 33 gates, the circuit of the 



square graph consists of only 18 gates, reflecting the simple symmetry of this graph. 



5.4. Random walk on the 8- star graph 



Fig. 10 shows another simple graph, namely the 8-star graph 
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Figure 10: The 8-star graph 



For this graph, the quantum walk operator is 



U =TC 



-0.75 0.25 0.25 0.25 0.25 0.25 0.25 0.25 \ 

0.25 -0.75 0.25 0.25 0.25 0.25 0.25 0.25 

0.25 0.25 -0.75 0.25 0.25 0.25 0.25 0.25 

0.25 0.25 0.25 -0.75 0.25 0.25 0.25 0.25 

0.25 0.25 0.25 0.25 -0.75 0.25 0.25 0.25 

0.25 0.25 0.25 0.25 0.25 -0.75 0.25 0.25 

0.25 0.25 0.25 0.25 0.25 0.25 -0.75 0.25 

0.25 0.25 0.25 0.25 0.25 0.25 0.25 -0.75 



which is a 16 x 16 unitary matrix. The complete result is displayed below. Note that since all of the rotation 
angles of A\, A\, A A V A\, A*, A 3 6 , A* are zeros, these gates vanish and are omitted here. 



Output of the package 



Matrix 



Gate 



GATEPI 

l; 

N 

GATEPI 

2;1 

NN 

GATEPI 
3; 1,2 
NNNN 

GATEPI 

4; 1,2, 3 

NNNNNNNY 



U 



16 C 3 n(4;(l,2,3) m ) 



- 7T - 
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Output of the package 



Matrix 



Gate 



Note: This part specifies the subset of gates: 




GATEY 
4; 1,2, 3 

0.0000 0.0000 0.0000 0.0000 
-0.1153 0.1573 -0.3073 -0.0372 



A 4 



C%(4; 1,2,3) 
29 5 = 2n- (-0.1153) 
26> 6 = 2tt- 0.1573 
26> 7 = 2n • (-0.3073) 
20 8 = In • (-0.0372) 



3 



Rl 



GATEY 

3; 1,2,4 

0.0000 0.0000 0.0000 0.0000 
-0.1289 -0.4029 0.0714 0.3210 



A 3 



C%(3; 1,2 A) 

26 5 = 2n -(-0.1289) 

29 6 = 2n • (-0.4029) 

29 7 = 2?r- 0.0714 
26> 8 = 2n ■ 0.3210 



3 



Rl 



GATEY 

4; 1,2, 3 

0.0000 0.0000 0.0000 0.0000 
0.0392 -0.2637 -0.4685 0.1927 



A 4 



C%(4; 1,2,3) 
29 5 = 2n ■ 0.0392 
26> 6 = 2n • (-0.2637) 
29 7 = 2n • (-0.4685) 
29* = 2tt- 0.1927 



-J. 



-Rl -Rt - Rl - Rl - 



GATEY 

2; 1,3,4 

0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 0.5000 



A 2 



C%{2; 1,3,4) 
29* = 2n-0.5 = n 



-R s , - 
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Output of the package 



Matrix 



Gate 



GATEY 
4; 1,2, 3 

0.0000 0.0000 0.0000 0.0000 
-0.0392 0.2637 0.4685 0.1927 



A 4 



C%(4; 1,2,3) 

20 5 = 2n • (-0.0392) 

20 6 = 2n • 0.2637 

20 7 = 2n ■ 0.4685 

20 8 = 2;r- 0.1927 



-J. 



-Rl -Rt -Rl - Rl - 



GATEY 
3; 1,2, 4 

0.0000 0.0000 0.0000 0.0000 
0.1289 0.4029 -0.0714 0.3210 



C%(3; 1,2,4) 
- 3 20 5 = 2tt- 0.1289 
ll ° 20 6 = 2n • 0.4029 

20 7 = 2n -(-0.0714) 
20« = 2n ■ 0.3210 



R 5 



tf 6 



GATEY 
4; 1,2, 3 

0.0000 0.0000 0.0000 0.0000 
0.1153 -0.1573 0.3073 -0.0372 



A 4 



C%(4; 1,2,3) 

20 5 = 2tt- 0.1153 

20 6 = 2tt ■ (-0.1573) 

20 7 = 2tt • 0.3073 

20 8 = 2n ■ (-0.0372) 



— 1 


I ( 


1 1 


» 1 

i 1 


p — 
i — 




) C 


) 1 




> — 


-R 5 V -R 

y 


T 

* -Rl -R 

y y 


8 

y 



GATEY 

1;2,3,4 

-0.5000 -0.5000 -0.5000 -0.5000 
-0.5000 -0.5000 -0.5000 -0.5000 



A 1 



C 3 i?,(l;2,3,4) 
20j = 202 = 203 = 

205 = 206 = 207 = 

2n ■ (-0.5) = -7T 



20 4 
20« 



RyU 



- Rl - Rl - 



-Rl - 



Putting these gates together, we obtain its quantum circuit as shown in Fig. [TT 



It requires only 34 quantum gates to implement the 16x16 unitary matrix. As a comparison, for a random 

16 x 16 real matrix, it needs around 8x15 + 15 = 135 gates, while for a random 16 x 16 complex matrix, it 
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Figure 11: The circuit of U, the random walk evolution on the star graph (Fig. 10 1 



needs 8x15 + 8x15 + 16 = 256 gates. The last 7 R y gate combinations (56 R y subgates) completely disappear. 
This is another example demonstrating the efficiency of Qcompiler. 

5.5. Random walk on a complicated graph 



Qcompiler also works for very complicated graphs. Fig. 12 is an adjacent matrix of a 100-node graph with 
white dots standing for '1' and black dots for '0'. 




Figure 12: The adjacent matrix of a complicated graph. 



The quantum walk operator for this graph is a 4011 x 4011 unitary matrix U . We first expand the size of 
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the unitary matrix with an identity matrix, i.e. 



W = 



U 



(32) 



where U is the 401 1 x 401 1 unitary matrix, / is a 85 x 85 identity matrix, and W is a 2 12 x 2 12 unitary matrix. 
The resulting quantum circuit contains 12 qubits. 



6. Conclusions 

We have developed an efficient and versatile package, Qcompiler, which maps any unitary matrix U of 
arbitrary size into a quantum circuit with only one- and two-qubit logic gates by applying cosine-sine de- 
composition recursively. For real unitary matrices, Qcompiler provides a much simpler quantum circuit with 
only half number of elementary gates in comparison with earlier work by Mottonen etal[9] which deals with 
general unitary matrices. 

The quantum circuits produced by Qcompiler also reflects the symmetry of the systems under study. In 
particular, we examined the resulting quantum circuits corresponding to quantum walks on graphs with certain 
degree of symmetry [19]. In this case, many gates in such a circuit turned out to be an identity gate or a simple 
n gate, which can be easily eliminated or combined with other gates to further reduce the complexity of the 
final quantum circuit. For unitary matrices of certain structure or symmetry, it may be possible to include 
qutris and qudits to produce simpler and more efficient quantum circuits. This will be an interesting subject 
for further study. 
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